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We measure the top quark mass in ti events using up to 3.6 fb _1 of pp collisions at y/s = 2 TeV. 
New results are described for matrix element methods with single lepton and dilepton chan- 
nels, and template methods with all-jets and dilepton channels. Standardization of sys- 
tematic uncertainties has proceeded between the CDF and DO collaborations. A combined 
measurement of 173. f ± 0.6(stat) ± l.l(syst) GeV is achieved. Fits to all measured elec- 
troweak parameters gives a Higgs boson mass 95% CL upper limit of tuh < 163 GeV. 
We also present a new measurement of the W boson mass from DO. Three different tem- 
plate methods are performed. Accounting for their correlations, a combined measurement of 
M w = 80.401 ± 0.043(stat + syst) GeV is obtained. 



1 Introduction 



The study of the top quark and W boson masses continues to be important in electroweak 
physics. Radiative corrections to the calculation of the W boson mass, My/, include terms 
which are logarithmic in the Higgs mass, m#, and vary as the square of the top quark mass, 
nit- This allows a constraint on run from measurements of the other two parameters. Another 
source of interest in the top quark mass has been the peculiar unitary value of the corresponding 
Yukawa coupling, Y t . As a result, the study of both the W boson and top quark masses are 
major goals of the Fermilab Tevatron program. I present recent measurements of both quantities 
in this paper. 

Experimentally, the measurements are currently made exclusively at the Fermilab Tevatron 
by the DO and CDF experiments. These are large, multipurpose collider detectors capable of low 
background particle identification for electrons and muons, and good momentum measurements 
for leptons, jets and $t- Recorded data samples amount to 6 fb" 1 per experiment. The mea- 
surement of nit has been undergoing rapid improvements in precision over the last few years. Up 
to 3.6 fb _1 have been analyzed so far, which represents a ten-fold increase from data available 



just four years ago. Improved methods have permitted more precise measurements. Advances 
in the estimation of systematic effects have substantially reduced remaining uncertainties. 

2 Methods of Top Quark Mass Measurement 

In general, the methods used to measure the top quark mass fall into two categories (for excep- 
tion, see Ref.^). Template methods fit a set of quantities, Xi, derived from event observables 
such as the measured kinematic quantities. The Xi are correlated with mt and are compared 
to distributions expected from top quarks of varying mass. Often, template methods involve a 
kinematic reconstruction of the event by solving constraint equations, and then fitting to the 
solved top quark mass. Matrix element methods instead compare the observables, the Xi of this 
method, directly to expectations using a leading order (LO) matrix element calculation for ti. 
In this way, event-by-event likelihoods are generated which use all of the kinematic information 
in an event. In each of these approaches, the background is accounted for by modeling the Xj. 
A probability of consistency vs. mt is then achieved by accounting for ti and background: 

P(xi\m t ) = fto P Ptop(xi\m t ) + fbkgPbkg (xi) (1) 

where Pt op and Pbkg are signal and background probability densities based on X{ and a particular 
mt- For matrix element approaches, Pt op incorporates the ti matrix element. For a particular 
final state channel: 

Pto P (xi\m t ) = — 1 r x f y2dq 1 dq 2 f(qi)f{q 2 )-^^==-d^e-W(x h y) 

Pobs(qq ->tt -» channel|m t ) J ^ Witl ' IV 

(2) 

where M is the LO matrix element, and f{q) are the parton distribution functions. 

The precision of the world average measurement of mt has become dominated by systematic 
uncertainties. The calibration of jet energies particularly limits the achievable precision of the 
mt measurement. However, top quark events offer an in situ dijet resonance in channels where 
W — > qq' — > jj. As a result, measurements in these channels typically perform a simultaneous fit 
to the jet energy scale (through m,,) and to mt- This allows the top quark mass to be evaluated 
at the scale favored by the correct M^y, and the residual systematic uncertainty is reduced. 
Systematic uncertainties have been defined to facilitate measurement combinations. Dominant 
physics modeling uncertainties come from hadronization, underlying event, and the modeling 
of backgrounds. Major detector modeling uncertainties arise from the residual jet energy scale, 
and also the difference between 6-jet and light-quark jet calibrations. Template methods can 
have significant uncertainty due to the statistical samples available for template generation. 

3 All-hadronic Channel 

In a ti event, both quarks are expected to decay to a W boson and a b quark. An all-hadronic 
(6j) final state results when both W bosons decay to quark pairs. Such a final state has 
an enormous background from QCD multijet production. This is ameliorated with the use 
of multiparameter discriminants and the tagging of 6-jets by identification of their displaced 
decay vertices from long b hadron lifetimes. CDF has performed a mass analysis in 2.9 fb _1 of 
collider data using this channel 1 ^. To enhance their sensitivity, they have developed a further 
discriminant which separates quark jets from gluon jets. The resulting background reduction 
gives a 25% improvement in statistical uncertainty from the mass measurement. 

The method of mass measurement proceeds by calculating dijet and three jet masses for 
each of the different assignments of jets to W bosons and b quarks. A x 2 is calculated of these 
masses with respect to fitted top quark and W boson masses, as well as fitted momenta, given 
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Figure 1: Left: CDF reconstructed mt from 6j events. Fitted ti and background contributions are overlaid with 
data. Right: DO fitted jet energy scale vs. mt for ^+jets events. Contours of constant AC are shown. 

expected widths and resolutions of these parameters. The shapes of the reconstructed W boson 
and top quark masses are obtained by fitting shapes from the ti and background samples. These 
are used as templates to compare to data. The result of this fit is shown in Fig. [TJ The use of 
the M\y constraint reduces the systematic uncertainty, but the residual jet energy scale is still 
a large contributor, along with color reconnection and residual bias errors. 



4 Single Lepton Plus Jets Channels 

Final states which provide high statistics but also reasonable background levels are those where 
at least one W boson from top quark decay itself decays to a lepton-neutrino (£v, where I = e, /i) 
pair. These £+jets events have been updated recently by both Tevatron experiments using matrix 
element approaches. DO has selected events in 2.6 fb -1 of Run 2b data using a neural network b- 
tag selection^. For both the signal and background probabilities given in Eq. [TJ matrix elements 
are used. As shown in Fig. [H the fitted jet energy scale is 1.018 db 0.007, which is consistent 
with the jet energy scale uncertainty from 7+jet events. The error is also within the expectation 
of pseudoexperiment tests from the Monte Carlo. The dominant systematic uncertainty comes 
from the relative 6/light jet response uncertainty of 0.8 GeV. 

Using 3.2 fb -1 , CDF has selected events with an exclusive four-jet selection^. One of these 
jets must satisfy a fo-tag requirement. A neural network selection is performed based on 19 
kinematic quantities using a quasi-Monte Carlo method to obtain uniformity in the integration. 
Not only is a clear separation between signal and background achieved, but the NN output is 
independent of top quark mass for the signal. In order to improve the mass resolution, the 
likelihood which comes from the mass fit is used to reject background and poorly reconstructed 
ti events. This is done by cutting out low peak likelihood events. The dominant uncertainty 
comes from the jet energy calibration. 

5 Dilepton Channels 

The rarest events occur when both W bosons decay to Iv. Such dilepton (21) channels can usually 
attain low backgrounds without resort to 6-tagging. However, there is no W — ► jj resonance 



Table 1: Recent m t measurements by CDF and DO experiments. The DO 2£ entry combines the e/x matrix element 
and ee/ nn/i+tTack template measurements. The last two entries give the combined result for each experiment 

for Run 1 plus Run 2 data. 



Channel 


Experiment 


Luminosity 


m t (GeV) 


6j 


CDF 


2.9 fb^ 1 


174.8 ± 2.4(stat + JES^i^syst) 1 ^ 


^+jets 


DO 


3.6 fb- 1 


173.7 ± 0.8(stat) ± 1.6(syst)EI 


^+jets 


CDF 


3.2 fb" 1 


172.1 ± 0.9(stat) ± 0.7(JES) ± l.l(syst)Sl 


2£ 


DO 


1-3.6 fb" 1 


174.7 ± 2.9(stat) ± 2.4(syst) b d 


all 


CDF 


up to 3.2 fb" 1 


172.6 ± 0.9(stat) ± 1.2(syst)^ 


all 


DO 


up to 3.6 fb -1 


174.2 ± 0.9(stat) ± 1.5(syst)El 



with which to obtain an in situ calibration and reduce the jet energy scale uncertainty. In 
addition, there is an ambiguity of how the measured event ]frr corresponds to the two neutrino 
momenta. The dilepton analyses have been pursued in both template and matrix element 
approaches recently. Using efi events in 3.6 fb -1 of data, DO has employed matrix elements for 
both signal and background in the likelihood calculation^. For the latter, the largest Z — > tt 
background is assumed for all backgrounds. This yields a measurement of m t = 174±3.3(stat) ± 
2.6(syst) GeV. The main systematic uncertainties are the jet energy scale uncertainty, and 
the difference between 6-jet and light quark jet energy scales. Hadronization and background 
modeling also play significant roles. 

Template methods have also been recently finalized by DO for the three standard dilepton 
channels (ee, e//, fi/j,) plus two 6-tagged channels with loosened lepton identification requirements 
(e//i+track)^. Because the dilepton event is underconstrained due to the two neutrinos, the 
kinematic reconstruction supplies an input parameter. They sample assumed neutrino rapidity 
distributions to provide a relative weight vs. mj. DO uses the moments of the weight distribution 
in 1 fb -1 to extract mt = 176.2±4.8(stat)±2.1(syst) GeV. DO has also applied a template method 
with the same sample using partial production and decay information. These two measurements 
combined give m t = 174.8 ± 4.4(stat) ± 2.0(syst) GeV. 

6 World Average and Electroweak Fits 

The measured value of mt for 6j, ^+jets and 21 channels is given in Table [TJ incorporating 
combinations of Run 2a and 2b results. Combined measurements for each experiment using all 
Run 1 and Run 2 results are also indicated. The world average top mass determined from these 
is 173.1 ± 0.6(stat) ± l.l(syst) GeV^, which corresponds to Y t = 0.995 ± 0.007. In the context 
of the electroweak model, fitting all measured electroweak parameters including mt yields a new 
constraint of the Higgs boson mass of mu < 163 GeV^l^ at 95% CL. The most likely value is 
m H = 90^27 GeV (see Fig. ED, to be compared with the direct LEP2 lower limit of 114.4 GeV. 

7 Measurement of the W Boson Mass 

The DO experiment has completed a preliminary, precision measurement of the W boson mass. 
They use 499,830 W — > ev events obtained from 1 fb -1 of collider data to perform three template 
analyses in parallel. These use the electron pr, the $t an d the transverse mass calculated from 
these two, mj- = y / 2p|,pj , (l — cos((/> e — This measurement requires a precise calibration 

of the electromagnetic calorimeter, which is provided by the 18,725 Z — > ee events in this data 
sample. A detailed modeling of the detector response to the recoiling system is also needed. 




Figure 2: Left: Fit to van using electroweak constraints including the new world average rat measurement. The 
BIB < 163 GeV limit, as well as the direct 95% c.l. exclusion band of 160 — 170 GeV masses from Tevatron 
searches are indicated. Right: The fitted distribution of mr in W — > ev events in 1 fb _1 of DO data. Residuals 
of the data relative to the fitted W boson template are shown at lower right. 



The analysis proceeds by comparing the distribution of the template variables with high 
statistics distributions from models of signal and background. For signal, these are obtained 
in 10 MeV steps using RESBOsEQ an d PHOTOS^. These templates must have 10 8 events 
each, so a tuned fast simulation of the detector is used. The Z — > ee process is modeled in 
a similar fashion to facilitate a carryover of the electromagnetic calibration to the W events. 
A parametrized functional form for the Z boson reconstructed mass distribution is taken from 
GEANTed Monte Carlo samples and tuned to the data. The fit of the data to the templates 
in the mass range of 70 to 110 GeV yields the electron energy calibration scale and offset 
with x 2 /dof = 153/160. Z boson events are also used to study the hadronic recoil to the W 
boson. A GEANTed Monte Carlo sample again provides the functional form for the detector 
response to the hard component. Contributions due to the spectator quarks and additional pp 
collisions are modeled by fitting to Z boson collider data. The backgrounds in this channel arise 
from Z — > ee events where one electron is not identified, QCD instrumental background, and 
W — ► tv — > evv. These are modeled using an electron plus track control sample, a trigger sample 
with track requirements omitted, and a full GEANT Monte Carlo simulation, respectively. 

The W boson mass is measured by blinding the data sample with an unknown offset for 
all three measurements simultaneously. The blinding was done to avoid bias by knowledge of 
the current world average value. The comparison of the data to the templates is performed, 
yielding good x 2 f° r each measurement (e.g. uit yielded x 2 /dof = 48/49). Once the blinding 
was removed, a measurement was obtained for each method, as given in Table [2] and Fig. [2j 
The systematic uncertainties are dominated by a 34 MeV uncertainty due to the electron energy 
scale. The methods are somewhat decorrelated, and a combination yields 80.401 =L 0.021 (stat) ± 
0.038(syst) GeV = 80.401 ± 0.043 GeVES. This will be propagated into the world average W 
boson mass and electroweak fits soon. 



8 Conclusions 



Tevatron experiments have analyzed up to 3.6 fb of collider data for top quark and W boson 
mass measurements. For the former, the all-jets channels are achieving 1.5% precision, the 
single leptons have surpassed the 1% precision, and even the rare dilepton events are providing 



Table 2: Measurements of Mw using three different template variables. W — > ev events in 1 fb 1 of DO data 

were used. 



template variable 


M w (GeV) 


rriT 


80.401 ± 0.023(stat) ± 0.037(syst) 




80.400 ± 0.027(stat) ± 0.040(syst) 


Pt 


80.402 ± 0.023(stat) ± 0.044(syst) 



2% precision to the top mass measurement. The world average is now 173.1 ± 0.6(stat) ± 
l.l(syst) GeV, which gives 0.7% precision. This update generates a new 95% CL upper limit 
on the Higgs boson mass of 163 GeV. The W boson mass has also been measured by DO in 
the ev channel, yielding the world's best measurement by a single experiment: Mw = 80.401 ± 
0.043(stat + syst) GeV. 
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